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Abstract. 

The region around PSR B1706-44 has been observed with the H.E.S.S. imaging atmospheric Cherenkov telescopes 
in 2003. No evidence for 7-ray emission in the TeV range was found at the pulsar position or at the radio arc 
which corresponds to the supernova remnant G 343.1-2.3. The 99 % confidence level flux upper limit at the pulsar 
position is F u j(E>350 GeV) = 1.4 x 10~ 12 s -1 cm~ 2 assuming a power law (dN/dE cx E~ r ) with photon index of 
r=2.5 and F m (E>500GeV) = 1.3 x 10~ 12 s _1 cm -2 without an assumption on the spectral shape. The reported 
upper limits correspond to 8% of the flux from an earlier detection by the CANGAROO experiment. 
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1. Introduction 

PSRB1706-44 is a young pulsar (spin-down age of 
^17kyrs) with distance estimates ranging from 1.8 to 
3.2 kpc with a period of 102 ms and a spin-down lumi- 
nosity of about 1% of the Crab pulsar (3.4 x 10 36 ergs -1 ). 
Pulsed emission has been observed at radio and X-ray 
wavelengths, and in GeV 7-rays. An extended synchrotron 
nebula around this compact object has been found in radio 
observations (IjGiacani et al. 20f)2j) with an extension of 1'— 
4' and with a flat spectrum (energy index of 0.3), and also 
in X-rays HGotthelf et al. 2002(1 with an extension of ^20" 
and with a photon index of 1.34. These characteristics sug- 
gest the existence of a pulsar wind nebula (PWN) pow- 
ered by the pulsar. In the TeV range, the CANGAROO 
experiment detected a steady emission coincident with the 
PWN position at a level of roughly 50 % of the Crab flux 
(Kifune et al. 119951 IKushida et al. 2003J) . suggesting that 
this PWN is the southern equivalent of the Crab nebula. 
The Durham Mark 6 collaboration (Chadwick et al. 11998(1 
reported also a significant detection above 300 GeV. A flux 
upper limit above 500 GeV which is compatible with the 
CANGAROO flux has been derived using data from the 
BIGRAT telescope (Rowell et al.^MB- 

PSRB1706-44 is coincident with an incomplete 
arc of radio emission l|McAdam et al. 1993|l which has 
been interpreted as a shell-type supernova remnant 
(SNR) named G 343.1-2.3. This SNR has been de- 
tected only at radio wavelengths l(Duncan et al. 1 9951 
and may be associated with the pulsar as discussed in 
Bock & Gvaramadze (2002) 

We present here the results of the observation 
of the field of view around PSR B1706-44 with 
the H.E.S.S. experiment. H.E.S.S. is an atmospheric 
Cherenkov detector dedicated to the observation of 
TeV 7-rays (Hofman n 2003(1 . Situated in Namibia, the 
full four-telescope array is operational since December 
2003. Each telescope has a mirror area of 107 m 2 
( Bernloh r et al. 20 03 1 and is equipped with a cam- 
era consisting of 960 photomultiplier tubes (PMT) 
l| Vincent et al. 2003(1 . The system has a field of view of 5°. 
In stereoscopic observation mode, it allows one to recon- 
struct the direction of individual showers with a precision 
better than 0.1°. 

2. Observations and data analysis 

PSR B1706-44 was observed with two H.E.S.S. telescopes 
between April and July 2003. During this commission- 
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Fig. 1. 9 2 distribution calculated with respect to the PWN 
position. The dots denote events from the ON region, the 
histogram are the events from the OFF region scaled by 
the normalization factor a. The dashed vertical line indi- 
cates the applied angular cut. The inset shows the differ- 
ence between the ON and the scaled OFF regions. 

ing phase, GPS time stamps were used in the offline data 
analysis to identify showers observed in coincidence by the 
two telescopes. This coincidence requirement allows for a 
higher background rejection and thus for a better sensi- 
tivity than single telescope observations. In this configu- 
ration, a source with a flux of 5% of Crab nebula can be 
detected with more than 5 a in 4.5 hours at 20° zenith an- 
gle. The pulsar was observed with 28-min runs in wobble 
mode, whereby runs are taken pointing ±0.5° away from 
the pulsar position in declination. Data affected by hard- 
ware problems or bad weather conditions were excluded 
from analysis. The proper functioning of the detector sys- 
tem was verified by numerous checks. The telescope point- 
ing has been confirmed by correlating high PMT currents 
with bright stars in the field of view. The trigger rate of 
the system is well reproduced by simulations for cosmic 
rays, and the shape of simulated 7-ray images is consis- 
tent with the result of Crab observations. Data analysis is 
performed with two completely independent chains with 
different calibrations, with independent Monte Carlo sim- 
ulations and with different analysis techniques. 

The selected data have a total live time of 14.3 hours. 
The energy threshold estimated from Monte Carlo sim- 
ulations at the average observation zenith angle (~26°) 
is about 350 GeV. This threshold is higher than for 
the four-telescope system since the telescopes were op- 
erated with higher trigger thresholds in the commission- 
ing phase. Data were analysed using standard shower 
reconstruction and standard background rejection meth- 
ods l(Aharonian et al. 2004(1 . Standard cuts, optimised on 
Monte Carlo simulations, have been applied on mean 
scaled Hillas parameters in order to increase the signal- 
to-background ratio. Showers were classified using the an- 
gular distance 9 between their reconstructed direction and 
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the direction of possible source. For this standard analy- 
sis, showers were accepted as coming from the source (the 
ON region) when their 9 2 was smaller than 0.02 degree 2 
(i.e. angular distance smaller than 8.5'). The background 
was determined by counting events in a ring (the OFF 
region) centered at the investigated direction whose inner 
radius is larger (>0.4°) than the Q 2 cut and whose area is 
7 times larger than the ON region. A normalization fac- 
tor a is applied to these estimated background counts to 
correct for the different size of ON and OFF regions and 
the different radial acceptance in the field of view. 

3. Results 

A plot of 2 relative to the PWN position is shown 
in Figure ^ The significance, calculated according to 
Li k Ma (1983)1 is 0.1 cr. Table □ provides an overview of 
the event statistics in the column labelled Standard. The 
analysis described above was repeated with the same cuts 
for every point in the field of view. The resulting signifi- 
cance map is presented in Fig. El It exhibits no significant 
point source excess in the vicinity of the pulsar or on the 
radio emission arc. 

In order to roughly reproduce the conditions of the 
PSR B1706-44 detection by CANGAROO, the analysis 
at the pulsar position was repeated using a looser 9 2 
cut of 0.05 degree 2 and selecting events above an energy 
of 1 TeV. The results are shown in the column labelled 
CANGAROO of Tableland give no indication for a sig- 
nificant excess. For the analysis of the radio arc, a 9 2 
cut of 0.36 degree 2 has been applied around the position 
(17h08m,— 44°17') and no significant excess is measured 
(column labelled Radio arc of Tabic 
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Standard 


CANGAROO 


Radio arc 


N m 


352 


112 


4746 


N oB 


2243 


512 


13688 


a 


0.15620 


0.19258 


0.34592 


Excess 


1.6 ± 20.2 


13.4 ±11.1 


11.0 ±79.9 


Significance 


O.lcr 


1.2 cr 


O.lcr 



Table 1. Analysis results: N on and N g are the event 
numbers in the ON and OFF regions, a is the normal- 
isation factor. The results are reported for the standard 
2 cut (column labelled Standard), for the cuts reproduc- 
ing the conditions of the CANGAROO detection (column 
CANGAROO) and for the analysis of the radio arc (col- 
umn Radio arc). 

Limits on the integral flux above certain energies Et 
were obtained using two different methods. The first 
method (Method A) tests the hypothesis that the num- 
ber of excess events with energies above Et result from a 
source with a power law spectrum with a (positive) pho- 
ton index T. The photon index was varied between 2 and 
3. This range includes the value of T=2.5 from the earlier 
CANGAROO detection. The second method (Method B) 
makes no assumption about the source spectrum and cal- 
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Fig. 2. Significance map centered on PSR B1706-44. 
The cross marks the pulsar position. The contour 
lines correspond to the 2.2 GHz image of G 343-2.3 
IjDuncan et al. 1995J) . The solid circle indicates the inte- 
gration region of the Standard cuts, the dashed circle the 
CANGAROO cuts and the dot-dashed circle the Radio arc 
cuts. The significance distribution for the entire H.E.S.S. 
field of view is compatible with a Gaussian of mean —0.06 
and of sigma 1.09. 



culates the integrated flux F directly as the difference of 
the measured flux from the ON region and the flux of 
cosmic-rays from the OFF region: 

/N on N off \ 

\l— 1 2 — 1 / 

Here, T is the live time, and both sums on the ON 
and OFF regions run over all showers with reconstructed 
energies greater than Et- The effective areas (Ai) depend 
on the zenith angle and energy of each event, and a is 
the normalization factor. As Ai is determined using the 
reconstructed energy, the energy threshold should be in- 
creased such that the bias of the reconstructed energy is 
less than 10 %. The upper limits derived with both meth- 
ods were calculated using the unified approach of Feldman 
& Cousins 1)1998[1 and a confidence level of 99%. To com- 
pare the upper limits from Method B with a prediction, 
the investigated model spectrum must be integrated over 
all energies starting at Et- 

Table |21 gives the values of flux upper limits at 99 % 
confidence level for various cuts and methods; both meth- 
ods give similar results. With method A, the upper limit 
at the PWN position corresponds to ~1% of the flux from 
the Crab Nebula (at the same energy threshold) and the 
upper limit for the radio arc corresponds to ~5% of the 
flux from the Crab Nebula. The upper limit which repro- 
duces the experimental conditions of the CANGAROO 
experiment corresponds to ~8 % of the flux reported by 
that collaboration. 
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Energy (TeV) 

Fig. 3. Integral upper limits at 99 % CL for the flux from 
the PWN position (solid, dotted and dot-dashed line). The 
filled circle corresponds to the CANGAROO detection 
of Kifune et al. QMS) and the CANGAROO integrated 
flux (grey area) is calculated from the result of a broken 
power law fit to the 2000 and 2001 differential spectrum 
UKushida et al. 2003}) . The open diamond and the triangle 
are from Rowell et al. I|1998|l and Chadwick et al. 1|1998|) . 
respectively. 



Method A Method B 

Standard 1.4 x 10" 12 (0.35) 1.3 x 10" 12 (0.50) 
CANGAROO 6.4 x 10~ 13 (1.00) 7.7 x 10~ 13 (1.00) 
Radio arc 5.8 x 10~ 12 (0.35) 3.5 x 10~ 12 (0.50) 

Table 2. Flux upper limits at 99% confidence level for 
the pulsar position in s -1 cm -2 . The upper limits from 
Method A were calculated assuming a photon index of 
r=2.5. The numbers in parentheses are the energy thresh- 
olds (in TeV) for which the upper limits were determined. 

4. Discussion and Conclusions 

The reported upper limits on the flux of TeV 7-rays are 
roughly one order of magnitude lower than the reported 
CANGAROO flux and a factor of 55 lower than earlier 
limits (Rowell et al. fiMty . The CANGAROO observations 
were not contemporaneous with the H.E.S.S. observations, 
which raises the question of whether the TeV emission 
could be variable on a time scale of years. Such a vari- 
ability seems unlikely given our current understanding of 
PWN (Blon din et al. 2001(1 . Another potential reason for 
the discrepancy could be an object confusion along the 
line of sight. There are, however, no BL Lac objects or 
variable galactic TeV 7-ray emitters known around the 
pulsar. It has been pointed out (Aharonian et al. 119971 
IKushida et al. 2003)l that the high flux level reported by 
CANGAROO is surprising. Since the X-ray luminosity is 
about 0.01 % of that of the Crab PWN, the TeV 7-rays 
should be emitted from a much larger volume than the 
X-rays, according to the inverse Compton (IC) scenario. 

Using the H.E.S.S. flux upper limit above ITeV, a 
lower limit on the magnetic field can be derived from 



Eq. 6 of Aharonian et al. I|1997f) . This requires a mea- 
surement of the flux in the X-ray band from the same 
electron population that emits the hypothetical TeV ra- 
diation. Measurements by Chandra l|Gotthelf et al. 2002|) 
provide a flux from the PWN, excluding the point-like 
emission of the central source; however, their chosen anal- 
ysis region (radius less than 10") is smaller than the full 
extent of the PWN, for which Finley et al. (1998) found 
a best-fit exponential scale length of 27". The flux mea- 
sured by ASCA (Finley et al. I1998J) encompasses the en- 
tire PWN, but also includes the pulsar emission. To esti- 
mate the PWN flux, we used the ASCA spectrum but 
subtracted a point source contribution estimated from 
ROSAT HRI to be (43 ± 12)% (Finley et al. ITMSI) . yield- 
ing an unabsorbed flux of 5.5 x 10 _13 erg s" 1 cm~ 2 in 
the 2-10 keV band for the PWN. The lack of observable 
X-ray emission below about 0.5 keV due to interstellar 
absorption means that the electrons producing the ob- 
served X-rays have somewhat higher energy than those 
producing TeV 7-rays, and an extrapolation of the X-ray 
spectrum to lower energies is necessary. The spectral in- 
dex measured with ASCA, r=1.7+Q'4, is fully compati- 
ble with the more precise determination from BeppoSAX, 
T= 1.69 ±0.29 (|Mineo et al. 2002jl . The derived lower limit 
on the magnetic field strength is then about 1 /J.G when 
one assumes that the inverse Compton scattering involves 
only the photons of the microwave background radiation 
and assuming the same photon index in the X-ray and TeV 
band. This value is however not very constraining given 
that the mean Galactic magnetic field is of the same order 
of magnitude. 
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